NON-GAUSSIAN  PROPERTIES  OF  WAVES  IN  FINITE  WATER  DEPTH 


By 


DAVID  J.  ROBUXARD 


A  THESIS  PRESENTED  TO  THE  GRADUATE  SCHOOL 

OF  THE  UNIVERSITY  OF  FLORIDA  IN  PARTIAL  FULFILLMENT 

OF  THE  REQUIREMENTS  FOR  THE  DEGREE  OF 

MASTER  OF  ENGINEERING 


UNIVERSITY  OF  FLORIDA 
1996 


Ip0 

'.I 


ACKNOWLEDGMENTS 

David  J.  Robillard  wishes  to  express  sincere  thanks  to  his  wife  and  children  for 
their  understanding  and  support.  The  author  would  also  like  to  thank  the  members  of  the 
graduate  committee,  Dr.  Michel  K.  Ochi,  Dr.  Robert  G.  Dean  and  Dr.  Ashish  Mehta,  for 
their  input  and  guidance  during  this  study.  Special  thanks  are  extended  to  the  committee 
chairman,  Dr.  Ochi,  whose  instruction  and  insight  enabled  me  to  learn  a  great  deal  on  the 
topic  of  this  study. 


11 


TABLE  OF  CONTENTS 

ACKNOWLEDGMENTS ii 

LIST  OF  FIGURES v 

LIST  OF  TABLES viii 

ABSTRACT ix 

1  INTRODUCTION 1 

1.1    Purpose  of  Study 2 

2  LITERATURE  REVIEW 5 

3  BACKGROUND/STUDY  APPROACH 8 

3.1  Wave  Record  Data 8 

3.1.1  Data  Acquisition 8 

3.1.2  Wave  Profile  Transformations 9 

3.2  Nearshore  Zone  Characteristics 16 

3.2. 1  Discussion  of  Water  Depth  Regions  During  ARSLOE 16 

3.2.2  Discussion  of  the  Surf  Zone  During  ARSLOE 17 

3.3  Non-Gaussian  Probability  Density  Function 19 

3.3.1  Non-Gaussian  Probability  Density  Function  Background 19 

3.3.2  Probability  Density  Function  Application  to  Wave  Data 21 

3.3.3  Dimensional  Analysis  of  PDF  Parameters 22 

4  PROBABILITY  DENSITY  FUNCTION  ANALYSIS 24 

4.1  Probability  Density  Function  Applicability  Verification 24 

4.2  Gaussian/Non-Gaussian  Boundary  Definition 33 

5  TREND  ANALYSIS  OF  PARAMETERS 42 

5.1  Relationship  Between  a  and  Water  Depth 42 

5.2  Non-Dimensional  PDF  Parameter  Analysis 44 

5.2.1  Spatial  Analysis  of  Data 44 

5.2.2  Relationship  Between  aa,  ii*/g  and  g*/o 52 

5.2.3  Relationship  Between  aa  and  a/d 54 

6  SUMMARY  AND  CONCLUSIONS 56 


in 


APPENDIX 

COMPARISON  OF  ORIGINAL  AND  SMOOTHED  "a"  VALUES 58 

LIST  OF  REFERENCES 66 

BIOGRAPHICAL  SKETCH 68 


IV 


LIST  OF  FIGURES 

Figure  page 

1.1.     Wave  Profile  Examples:  a)  Gage  615,  Severe  Sea  Condition,  and  b)  Gage 
625,  Mild  Sea  Condition  1 

3.1.  Coastal  Engineering  Research  Center  Field  Research  Facility 10 

3.2.  Gage  615  -  10/23  @  2155  Wave  Profile 12 

3.3.  Gage  635  -  10/23  @  2155  Wave  Profile 12 

3.4.  Gage  645  -  10/23  @  2155  Wave  Profile 12 

3.5.  Gage  655  -  10/23  @  2155  Wave  Profile 12 

3.6.  Gage  665  -  10/23  @  2155  Wave  Profile 13 

3.7.  Gage  675  -  10/23  @  2155  Wave  Profile 13 

3.8.  Gage  625  -  10/23  @  2155  Wave  Profile 13 

3.9.  Gage  615  -  10/25  @  1315  Wave  Profile 14 

3.10.  Gage  635-  10/25  @  1315  Wave  Profile 14 

3.11.  Gage  645-  10/25  @  1315  Wave  Profile 14 

3.12.  Gage  655  -  10/25  @  1315  Wave  Profile 14 

3.13.  Gage  665  -  10/25  @  1315  Wave  Profile 15 

3.14.  Gage  675  -  10/25  @  1315  Wave  Profile 15 

3.15.  Gage  625  -  10/25  @  1315  Wave  Profile 15 

4.1.  Gage  615-  10/23  @  2155  Histogram  Data 25 

4.2.  Gage  635-  10/23  @  2155  Histogram  Data 26 

4.3.  Gage  645-  10/23  @  2155  Histogram  Data 26 


4.4.  Gage  655  -  10/23  @  2155  Histogram  Data  27 

4.5.  Gage  665  -  10/23  @  2155  Histogram  Data  27 

4.6.  Gage  675  -  10/23  @  2155  Histogram  Data  28 

4.7.  Gage  625-  10/23  @  2155  Histogram  Data  28 

4.8.  Gage  615-  10/25  @  1315  Histogram  Data  29 

4.9.  Gage  635  -  10/25  @  1315  Histogram  Data  29 

4.10.  Gage  645-  10/25  @  1315  Histogram  Data  30 

4.11.  Gage  655-  10/25  @  1315  Histogram  Data  30 

4.12.  Gage  665-  10/25  @  1315  Histogram  Data  31 

4.13.  Gage  675-  10/25  @  1315  Histogram  Data  31 

4.14.  Gage  625-  10/25  @  1315  Histogram  Data  32 

4.15.  Sample  Data  Points  to  Analyze  Gaussian/non-Gaussian  Boundary 34 

4. 16.  Probability  Density  Function  of  POINT  A  of  Figure  4. 15 35 

4.17.  Probability  Density  Function  of  POINT  B  of  Figure  4. 15 36 

4.18.  Probability  Density  Function  of  POINT  C  of  Figure  4. 15 36 

4. 19.  Probability  Density  Function  of  POINT  D  of  Figure  4. 15 37 

4.20.  Probability  Density  Function  of  POINT  E  of  Figure  4. 15 37 

4.21.  Probability  Density  Function  of  POINT  F  of  Figure  4. 15 38 

4.22.  Probability  Density  Function  of  POINT  G  of  Figure  4. 15 38 

4.23.  Probability  Density  Function  of  POINT  H  of  Figure  4. 15 39 

4.24.  Gaussian/N on- Gaussian  Boundary  Defined 41 

5.1.  Spatial  Comparison  of  a  at  Mild  and  Severe  Sea  Conditions 43 

5.2.  All  c  Data  Plotted  as  a  Function  of  Water  Depth 43 

5.3.  Spatial  Comparison  of  ao  at  Mild  and  Severe  Sea  Conditions 45 

5.4.  Spatial  Comparison  of  u*/a  at  Mild  and  Severe  Sea  Conditions 46 


VI 


5.5.  Spatial  Comparison  of  o*/a  at  Mild  and  Severe  Sea  Conditions 47 

5.6.  aa  as  a  Function  of  Water  Depth 49 

5.7.  u*/a  as  a  Function  of  Water  Depth 50 

5.8.  ct*/q  as  a  Function  of  Water  Depth 51 

5.9.  Plot  of  |Wa  as  a  Function  of  aa  for  Gages  615  and  625 53 

5.10.  Plot  of  g*/g  as  a  Function  of  ao  for  Gages  615  and  625 53 

5.11.  Plot  ofaa  as  a  Function  of  o/d  for  Mild  and  Severe  Sea  Conditions 55 

A.l.  Plot  of  Gage  615  Original  and  Smoothed  "a"  Values 60 

A.2.  Plot  of  Gage  635  Original  and  Smoothed  "a"  Values 61 

A.3.  Plot  of  Gage  645  Original  and  Smoothed  "a"  Values 62 

A.4.  Plot  of  Gage  655  Original  and  Smoothed  "a"  Values 63 

A.5.  Plot  of  Gage  665  Original  and  Smoothed  "a"  Values 64 

A.6.  Plot  of  Gage  675  Original  and  Smoothed  "a"  Values 65 

A.7.  Plot  of  Gage  625  Original  and  Smoothed  "a"  Values 66 


Vll 


LIST  OF  TABLES 
Table  page 

3.1.  Deep  Water  Wave  Statistics  for  Mild  and  Severe  Sea  Conditions 16 

3.2.  Surf  Zone  Locations  for  Mild  and  Severe  Sea  Conditions 19 

4.1.  Wave  Profile  Statistics  for  Gage  615  -  10/23  @  2155 25 

4.2.  Wave  Profile  Statistics  for  Gage  635  -  10/23  @  2155 26 

4.3.  Wave  Profile  Statistics  for  Gage  645  -  10/23  @  2155 26 

4.4.  Wave  Profile  Statistics  for  Gage  655  -  10/23  @  2155 27 

4.5.  Wave  Profile  Statistics  for  Gage  665  -  10/23  @  2155 27 

4.6.  Wave  Profile  Statistics  for  Gage  675  -  10/23  @  2155 28 

4.7.  Wave  Profile  Statistics  for  Gage  625  -  10/23  @  2155 28 

4.8.  Wave  Profile  Statistics  for  Gage  615-  10/25  @  1315 29 

4.9.  Wave  Profile  Statistics  for  Gage  635  -  10/25  @  1315 29 

4.10.  Wave  Profile  Statistics  for  Gage  645  -  10/25  @  1315 30 

4.11.  Wave  Profile  Statistics  for  Gage  655  -  10/25  @  1315 30 

4.12.  Wave  Profile  Statistics  for  Gage  665  -  10/25  @  1315 31 

4.13.  Wave  Profile  Statistics  for  Gage  675  -  10/25  @  1315 31 

4.14.  Wave  Profile  Statistics  for  Gage  625  -  10/25  @  1315 32 

4.15.  Statistical  Properties  of  Data  Points  A-H  in  Figure  4.15 35 

5. 1.     General  Shallow  Water  Limits  of  Parameters 52 


Vlll 


Abstract  of  Thesis  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Master  of  Engineering 


NON-GAUSSIAN  PROPERTIES  OF  WAVES  IN  FINITE  WATER  DEPTH 

By 

David  J.  Robillard 

May  1996 

Chairman:  Michel  K.  Ochi 

Major  Department:  Coastal  and  Oceanographic  Engineering 

As  waves  propagate  into  shallower  waters,  their  profile  changes  from  a  Gaussian 
distribution  to  an  increasingly  non-Gaussian  profile.  The  non-Gaussian  properties  of 
waves  in  finite  water  depths  are  analyzed  through  the  use  of  a  probability  density  function 
defined  in  closed  form  and  that  can  be  solved  by  the  application  of  wave  displacement 
data.   Initially,  the  applicability  of  the  probability  density  function  is  verified  using  wave 
data  records  obtained  from  the  Coastal  Engineering  Research  Center  (CERC)  at  Duck, 
North  Carolina.  A  broad  range  of  sea  conditions  and  water  depths  are  represented  to 
ensure  the  probability  density  function  correlates  well  with  actual  histogram  data  over 
varied  conditions.  With  the  applicability  verified,  the  probability  density  function  is  used 
to  define  criteria  for  the  boundary  where  the  wave  field  can  no  longer  be  considered 
Gaussian.  The  ability  to  determine  when  a  shoaling  wave  profile  can  no  longer  be 
considered  Gaussian  is  of  paramount  importance  in  order  to  ensure  proper  wave  theories 
are  being  applied.  This  boundary  between  Gaussian  and  non-Gaussian  wave  profiles  will 
be  determined  in  terms  of  significant  wave  height  and  water  depth. 


IX 


Presently,  wave  data  records  are  required  to  define  the  probability  density  function 
used  in  this  study.  Trend  analyses  were  conducted  on  the  parameters  of  the  probability 
density  function  to  determine  if  simplifications  can  be  made  to  its  form  and  to  see  if  it  is 
possible  to  define  the  probability  density  function  in  terms  of  the  local  conditions,  sea 
severity  and  water  depth,  without  the  need  for  wave  data  records. 


CHAPTER  1 
INTRODUCTION 

As  wind  generated  waves  transition  from  deep  water  to  shallow  water,  a  significant 

change  is  evident  in  the  wave  profile.  As  an  example,  Figure  1-1  shows  portions  of  wave 

records  obtained  by  the  Coastal  Engineering  Research  Center,  (CERC),  during  the 

ARSLOE  Project.  Figure  l-l(a)  is  the  wave  profile  recorded  at  a  distance  of  60  meters 

offshore  in  a  shallow  intermediate  water  depth  condition.  Figure  1-1  (b)  is  the  wave 

profile  taken  at  the  same  time  but  at  a  distance  of  12  km  offshore  in  a  deep  water 

condition.  Deep  water  wave  profiles  in  the  sea  characteristically  reflect  a  Gaussian  (linear) 

random  process.  The  height  of  the  wave  crests  are  probabilistically  equal  to  the  height  of 

the  wave  troughs.  As  waves  move  from  deep  to  shallow  water,  the  wave  profile 

transforms  from  a  Gaussian  random  process  to  a  non-Gaussian  (nonlinear)  random 

process  in  which  the  wave  profile  shows  a  definite  excess  of  high  crests  and  shallow 

troughs. 
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a)  Gage  615,  Severe  Sea  Condition  b)  Gage  625,  Mild  Sea  Condition 


Figure  1.1  Wave  Profile  Examples 


This  remarkable  transformation  in  the  wave  profile,  as  an  irregular  sea  transitions 
from  deep  to  shallow  water,  has  also  been  evidenced  in  laboratory  experiments  (Doering 
and  Donelan  1993).  The  transformation  of  wave  profiles  can  be  attributed  to  the  non- 
linear wave-wave  interactions  augmented  primarily  by  decreasing  water  depth.  The 
"degree"  of  non-linearity  a  wave  field  has  at  a  particular  site  of  interest  in  the  coastal 
region  appears  to  depend  largely  on  sea  severity  and  water  depth. 

Although  there  has  been  awareness  that  Gaussian-based  solutions  do  not  always 
represent  accurate  results  in  the  coastal  region,  clarification  of  wave  characteristics  from 
deep  to  shallow  water  has  been  a  subject  of  considerable  interest  in  ocean  and  coastal 
engineering.  There  have  been  few  studies  which  have  offered  non-Gaussian  solutions 
describing  the  probability  density  function  of  waves  in  coastal  regions.  A  commonly 
known  study  was  by  Longuet-Higgins  (1963),  which  statistically  predicts  coastal  wave 
properties  using  a  power  series  expansion.  This  probability  distribution  is  based  essentially 
on  the  Gram-Charlier  series  distribution.  A  probability  density  function  developed  by  Ochi 
and  Ann  (1994a)  appears  to  be  applicable  throughout  the  coastal  region  from  deep  to 


shallow  water  and  is  defined  in  closed  form  Equation  (1) 
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f(x)  = — exp 
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This  is  the  probability  density  function  that  will  be  used  as  the  basis  to  analyze 
coastal  wave  profiles  for  this  study. 

1.1    Purpose  of  the  Study 

The  overall  objective  of  this  study  involves  two  components.  The  first  component 
will  be  to  validate  the  applicability  of  the  probability  density  function  developed  for  waves 
in  finite  water  depth  (Ochi  and  Ahn  1994a)  and  to  analyze  the  transition  point  where  wave 


profiles  can  no  longer  be  considered  Gaussian.  In  order  to  accomplish  the  first  objective  of 
this  study,  a  statistical  analysis  was  carried  out  on  wave  time  histories  obtained  during  a 
storm  by  the  Coastal  Engineering  Research  Center  (CERC)  during  the  Atlantic  Ocean 
Remote  Sensing  Land-Ocean  Experiment  (ARSLOE)  Project  in  1980.  To  interpret  the 
transformation  of  waves  from  a  Gaussian  to  a  non-Gaussian  random  process,  wave 
records  at  seven  locations  were  taken  where  the  water  depths  ranged  from  two  meters  to 
twenty-one  meters.  In  total,  approximately  1000  wave  records  were  analyzed  with  each 
wave  record  consisting  of  4800  water  displacement  measurements.  The  individual  wave 
displacement  measurements  were  recorded  at  a  frequency  of  4  Hertz  in  twenty  minute 
intervals. 

Upon  confirmation  of  the  probability  density  functions'  applicability,  the  histogram 
data  will  be  compared  closely  with  the  corresponding  normal  distribution  to  determine  the 
conditions  under  which  the  wave  profile  can  no  longer  be  considered  Gaussian.  In 
particular,  general  relationships  to  sea  severity  and  water  depth  will  be  determined.  A 
criteria  will  be  established  which  will  define  the  limit  to  when  the  Gaussian  wave  profile 
becomes  non-Gaussian.  This  limit  will  be  based  on  the  nonlinearity  parameters  analyzed 
from  the  wave  profile  data.  The  second  objective  of  the  study  will  then  be  to  analyze  the 
parameters  of  the  probability  density  function.  The  three  parameters,  a,  u*  and  o*  will  be 
statistically  analyzed  in  both  dimensional  and  non-dimensional  form.  Trends  based  on  the 
analysis  will  be  described  and  relationships  linking  these  parameters  of  the  probability 
density  function  to  site  conditions  will  be  developed. 

This  study  consists  of  six  chapters.  Chapter  2  provides  a  general  review  on  the 
development  of  probability  density  functions  applicable  to  coastal  regions.  Chapter  3 
describes  the  ARSLOE  Project  of  1980  and  addresses  the  statistical  approach  taken  in  this 
study.  The  data  used  as  the  basis  for  this  study,  were  obtained  from  the  ARSLOE  Project. 
This  chapter  will  also  describe  the  local  field  conditions  at  CERC  during  the  time  the  data 
were  collected.  Chapter  4  addresses  the  probability  density  function  applicability  and 


analysis.  This  chapter  has  two  distinct  sections.  The  first  section  reviews  the  finite  water 
depth  probability  density  function  (Ochi  and  Ann  1 994b)  and  verifies  the  applicability  of 
the  probability  density  function.  The  second  section  discusses  sea  severity  and  water 
depth  and  the  transition  between  Gaussian  and  non-Gaussian  wave  profiles.  Chapter  5 
relates  the  non-dimensionalized  parameters  of  the  probability  density  function,  aa,  u»/a 
and  (5*1(5,  to  sea  severity,  (a),  and  water  depth  in  both  dimensional  and  non-dimensional 
form.  The  parameters  will  also  be  compared  to  each  other  and  relationships  determined. 
This  chapter  describes  trends  evident  based  on  statistical  analysis.  These  trends  are 
qualitative  and  are  intended  to  give  a  better  understanding  of  the  interaction  the  wave 
profile  has  with  its  environment  as  it  propagates  in  shallow  waters.  Finally,  Chapter  6 
provides  the  summary  and  conclusion  for  this  study. 


CHAPTER  2 
LITERATURE  REVIEW 


As  stated  in  the  Introduction,  the  profile  of  waves  in  finite  water  depths  is  quite 
different  from  that  observed  in  deep  water.  From  the  stochastic  view  point,  waves  in  finite 
water  depth  are  categorized  as  a  non-Gaussian  random  process  for  which  a  precise 
description  of  the  process  does  not  exist.  This  is  in  contrast  to  the  Gaussian  random 
process  which  can  be  defined  by  a  rigorous  mathematical  formulation.  Although 
considerable  attention  has  been  given  to  Gaussian  waves  in  deep  water,  little  information 
is  available  on  the  non-Gaussian  characteristics  of  shallow  water  waves. 

One  of  the  most  commonly  known  presentations  describing  non-Gaussian  waves  is 
the  Gram-Charlier  series  probability  density  function.  This  density  function  is  given  in  a 
series  consisting  of  an  infinite  number  of  terms.  Longuet-Higgins  (1963)  theoretically 
showed  that  the  statistical  properties  of  nonlinear  waves  can  be  presented  by  the  Gram- 
Charlier  series  distribution.  Bitner  (1980)  found  reasonable  agreement  in  a  comparison 
between  the  Gram-Charlier  series  distribution  and  the  histogram  constructed  from  wave 
data  obtained  in  a  relatively  shallow  water  area. 

On  the  other  hand,  Ochi  and  Wang  (1984)  analyzed  more  than  500  wave  records 
obtained  by  the  Coastal  Engineering  Research  Center  (CERC)  at  Duck,  North  Carolina 
during  the  ARSLOE  Project  (1980).  They  found  that  while  the  statistical  properties  of 
waves  at  any  water  depth  can  be  well  represented  by  the  Gram-Charlier  series  distribution, 
there  is  a  drawback  to  applying  this  density  function  in  practice.  The  probability  density 
function  becomes  negative  for  large  negative  wave  displacements  and  this  negative  portion 
does  not  vanish  with  increasing  number  of  terms  in  the  series  distribution.  Because  of  this 


difficulty,  the  transition  property  of  wave  characteristics  from  gaussian  to  non-gaussian 
distribution  can  not  be  well  clarified.  They  derived,  however,  a  limiting  sea  severity  for  a 
specified  water  depth  below  which  wind-generated  coastal  waves  may  be  considered  to  be 
a  Gaussian  random  process. 

Another  approach  for  describing  the  surface  elevation  of  non- Gaussian  random 
waves  is  to  assume  that  individual  waves  can  be  expressed  as  a  Stokes  wave.  There 
appears  to  be  some  reservation  to  applying  this  approach  for  representing  the  statistical 
distribution  of  the  wave  profile.  This  method  presupposes  a  preliminary  form  of  the  wave 
profile  before  establishing  its  probability  density.  This  supposition  may  impose 
unfavorable  bias  on  the  outcome. 

Tayfun  (1980)  derived  the  probability  density  function  based  on  the  Stokes  profile 
approximation  up  to  the  second  order.  Huang  et  al.  (1983)  developed  the  probability 
density  function  considering  all  terms  up  to  the  third  order.  These  probability  density 
functions  are  superior  to  the  Gram-Charlier  series  distribution  in  that  they  do  not  have  any 
negative  portion  over  the  entire  range  of  wave  displacement.  Comparison  of  these 
probability  density  functions  and  wave  data  show  good  agreement. 

Langley  (1987)  applied  a  probability  density  function  representing  the  response  of 
a  non-linear  system  developed  by  Kac  and  Siegert  (1947),  to  nonlinear  random  waves. 
This  probability  density  function  is  given  in  the  form  of  a  series  and  can  only  be  evaluated 
numerically.  Ochi  and  Ann  (1994a)  modified  the  Kac  Siegert  density  function  so  that  it  is 
presented  in  closed  form  with  parameters  determined  through  spectral  analysis.  They 
applied  the  probability  density  function  to  non-Gaussian  coastal  waves.   In  the  application 
process,  they  decomposed  the  spectra  into  linear  and  nonlinear  components  to  clarify  the 
degree  of  non-linearity  involved  in  shallow  water  waves  (Ochi  and  Ann  1994b).  The 
probability  density  function  does  not  have  any  negative  portion  in  its  entire  domain  and 
appears  to  represent  the  statistical  distribution  of  the  wave  profile  in  any  sea  state  during  a 
storm  in  deep,  intermediate  and  shallow  water. 


Studies  concerned  with  the  non- Gaussian  properties  of  waves  in  finite  water  depth 
include  Thornton  and  Guza  (1983)  who  claimed  that  the  results  of  their  analysis  of  field 
data  obtained  at  Torrey  Pines  Beach,  California,  show  that  the  Rayleigh  probability 
distribution  can  be  applied  to  represent  wave  heights  observed  from  shallow  to  deep  (10.8 
meter)  water.  This  implies  that  all  observed  waves  can  be  considered  a  Gaussian  random 
process.  However,  the  sea  severity  in  relation  to  water  depth  at  all  locations  in  the 
nearshore  zone  where  their  measurements  were  taken  was  too  mild  and  below  the  limit  for 
waves  to  be  a  non-Gaussian  random  process.  This  subject  will  be  discussed  in  detail  in 
this  present  study. 

In  closing  this  chapter,  it  is  important  to  note  that  there  still  remains  a  significant 
amount  of  work  to  do  in  refining  our  ability  to  accurately  describe  irregular  seas  in  finite 
water  depths  and  over  a  wide  range  of  sea  severity. 


CHAPTER  3 
BACKGROUND  /  STUDY  APPROACH 


3.1    Wave  Record  Data 


An  experiment  jointly  hosted  by  the  U.S.  Army  Coastal  Engineering  Research 
Center  (CERC)  and  the  National  Ocean  Survey  (NOS)  was  conducted  from  06  October 
1980  to  30  November  1980.  An  aspect  of  this  Atlantic  Ocean  Remote  Sensing  Land- 
Ocean  Experiment  (ARSLOE)  involved  ocean  wave  studies  implemented  near  the 
CERC's  Field  Research  Facility  (FRF)  at  Duck,  North  Carolina.  In  situ  wave  measuring 
devices  were  deployed  in  mean  water  depths  ranging  from  1.35  meters  to  24.4  meters  and 
were  placed  in  a  linear  array  extending  from  60  meters  offshore  to  12,000  meters  offshore, 
Figure  3.1. 

3.1.1    Data  Acquisition 

The  wave  displacement  data  analyzed  in  this  study  came  from  a  total  of  eight  wave 
measuring  devices.  The  FRF  has  a  550  meter  research  pier  oriented  perpendicular  to  the 
local  shore  which  was  equipped  with  seven  Baylor  resistance-type  wave  staffs.  These 
Baylor  wave  staffs  function  by  detennining  the  change  in  electrical  resistance  experienced 
by  changes  in  water  elevations.  Round-trip  travel  time  of  an  electrical  current  from  an 
above-water  emitter  to  the  water  surface  is  measured.    Following  along  the  line  of  the 
pier  off-shore  was  a  Waverider  Buoy  Gauge.  This  is  a  surface-following  class  buoy.  It 
functions  by  double-integrating  the  output  of  an  accelerometer  held  vertical  by  a  pendulum 
to  produce  the  time  history  of  buoy  heave  motion.  The  statistics  of  the  buoy  heave  then 
represent  the  wave  displacements.  In  the  Waverider  buoy,  the  pendulum  is  contained 


entirely  within  the  spherical  float  so  the  buoy  does  not  have  to  be  stabilized.  These 
devices  recorded  both  wave  displacements  and  wave  periods  but  only  the  displacement 
data  were  analyzed  for  this  study.  The  wave  gauges  recorded  water  surface  fluctuations 
simultaneously  at  all  eight  locations  at  a  frequency  of  4  Hertz  over  twenty  minute 
intervals.  This  resulted  in  4800  individual  wave  displacements  recorded  for  each  twenty 
minute  segments  continuously.  Due  to  the  significant  amount  of  data  obtained  during  this 
unique  experiment,  only  a  sub-set  of  the  data  was  analyzed.  This  chosen  data  sub-set 
began  at  1200  on  23  October  1980  and  ended  at  1315  on  25  October  1980.  This  two  day 
segment  was  unique  in  that  a  severe  storm  occurred  from  23-25  October  which  was  the 
most  significant  meteorological  event  recorded  during  the  experiment.  Wind  speeds  on 
the  order  of  fifteen  meters  per  second  were  recorded  during  the  severe  storm  This  two 
day  sub-set  of  data  was  used  as  the  basis  for  this  study  since  it  included  the  entire  range  of 
wave  profiles  experienced  at  each  of  the  eight  wave  gauge  locations.  From  this  two  day 
sub-set,  a  total  of  nearly  5  million  individual  wave  displacements  were  analyzed  and  the 
wave  profile  data  covered  the  entire  range  of  sea  conditions  present  during  the  complete 
56  day  ARSLOE  Project. 

3.1.2    Wave  Profile  Transformations 

The  data  collected  at  the  FRF  during  ARSLOE  is  unique  and  impressive.  It  is  very 
well  suited  for  analyzing  changes  in  the  wave  field  as  the  waves  propagate  through 
shoaling  depths.  The  bathymetry  at  the  experiment  location  is  relatively  simple  showing 
straight  and  parallel  contours.  This  provided  favorable  conditions  which  maximized 
control  over  the  environmental  impact  on  the  quality  of  the  data  received.  The  wave 
profile  data  provide  a  matrix  of  wave  profiles  allowing  both  spatial  and  temporal  analyses 
to  be  performed. 
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Figure  3.1  Coastal  Engineering  Research  Center  Field  Research  Facility 
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A  wide  range  of  wave  profiles  was  observed  during  the  experiment.  Figures  3.2- 
3.8  represent  a  mild  sea  condition  recorded  on  10/23  @  2155  for  all  finite  water  depth 
gage  locations  shown  in  Figure  3. 1.  Figures  3.9-3. 15  represent  a  severe  sea  condition 
recorded  at  the  height  of  a  storm  on  10/25  @  1315  for  the  same  gage  locations.  It  is 
evident  from  the  sample  wave  data  that  the  wave  field  transitions  as  it  shoals  from  deep  to 
shallow  water.  The  wave  data  representing  the  mild  sea  condition  at  Gage  625,  Figure 
3.8,  are  characterized  by  an  almost  symmetrical  distribution  of  positive  and  negative  wave 
profile  displacements.  This  is  characteristic  of  Gaussian  random  processes.  As  the  wave 
shoals  into  shallower  water,  the  wave  profile  transforms  showing  progressively  higher 
peaks  and  shallower  troughs.  This  process  is  clearly  non-Gaussian  and  must  be  analyzed 
differently.  It  is  also  clear  from  observing  these  wave  records  that  the  degree  of  non- 
linearity  increases  as  the  relative  water  depth  decreases  and  as  sea  severity  increases.  This 
implies  that  analysis  of  this  nonlinearity  must  include  a  dependence  on  water  depth  as  it 
relates  to  the  magnitude  of  the  waves'  displacement  (sea  severity). 
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Figure  3-2  Gage  615  -  10/23  @  2155  (d=1.74m) 
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Figure  3-3  Gage  635  -  10/23  @  2155  (d=2.83m) 
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Figure  3-4  Gage  645  -  10/23  @  2155  (d=4.31m) 
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Figure  3-5  Gage  655  -  10/23  @  2155  (d=6.0m) 
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Figure  3-6  Gage  665  -  10/23  @  2155  (d=6.08m) 


Figure  3-7  Gage  675  -  10/23  @  2155  (d=6.82m) 
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Figure  3-8  Gage  625  -  10/23  @  2155  (d=9.66m) 


14 


IjJJ^J  JvUiJV^J  J>  ^jMK^J  a, 


Figure  3-9  Gage  615  -  10/25  @  1315  (d=2.23m) 
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Figure  3-10  Gage  635  -  10/25  @  1315  (d=3.32m) 
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Figure  3-11  Gage  645  -  10/25  @  1315  (d=4.66m) 
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Figure  3-12  Gage  655  -  10/25  @  1315  (d=6.44m) 
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Figure  3-13  Gage  665  -  10/25  @  1315  (d=6.52m) 


<y* 


Figure  3-14  Gage  675  -  10/25  @  1315  (d=7.22m) 
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Figure  3-15  Gage  625  -  10/25  @  1315  (d=10.04m) 
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3.2    Nearshore  Zone  Characteristics 

All  wave  fields  propagating  in  finite  water  depths  are  influenced  by  the  attributes 
of  the  local  sea  environment.  These  attributes  include  water  depth  and  beach  slope.  An 
awareness  of  the  characteristics  of  these  attributes  is  important  to  understand  the  type  of 
influence  they  have  on  the  wave  field. 

3.2.1    Discussion  of  Water  Depth  Regions  During  ARSLOE 

Based  on  linear  theory,  there  are  three  water  depth  regions  which  are  defined  using 

the  velocity  potential  for  a  progressive  wave.  The  velocity  potential  is  defined  as 

x     HgTcosk(k(d  +  z))  .  n       2**r 

Kx>z>')  =        a un^        sm(** — )  • 

4;rcosh(Az/)  T 

The  distinct  water  depth  regions  are  determined  by  the  asymptotes  of  the  hyperbolic 
function.  Deep  water  wave  data  obtained  at  gage  location  710  and  shown  in  Table  3. 1 
allow  these  regions  to  be  defined  for  the  mild  sea  condition  and  the  severe  sea  condition. 
In  general,  the  deep  water  wave  condition  is  defined  as  kd  >  k,  where  k  is  the  wave 
number  and  h  is  the  water  depth.  This  condition  occurs  at  d  =  21m  for  the  mild  sea 
condition  on  10/23  @  2155  and  d  =  43m  for  the  severe  sea  condition  on  10/25  @  1315. 
The  shallow  water  condition  is  defined  as  kd  <  rc/lO.  This  condition  occurs  at  d  =  0.6m  for 
the  mild  sea  condition  and  d  =  1.3m  for  the  severe  sea  condition  on  25  October.  The 
intermediate  water  condition  occurs  between  these  bounds.    All  wave  profile  data 
obtained  from  ARSLOE  are  located  in  the  intermediate  water  depth  region. 

Table  3.1  Deep  Water  Wave  Statistics  for  Mild  and  Severe  Sea  Conditions 


Date/Time 


10/23  @  2155 
10/25  (gj  1315 


Zero  Crossing  Period,  T 


5.17  sec. 
7.43  sec. 


Deep  Water  Wavelength,  X.0 


41.715  m 
86.157  m 
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3.2.2    Discussion  of  the  Surf  Zone  During  ARSLOE 

An  integral  aspect  of  any  study  which  covers  the  nearshore  region  must  analyze 
and  consider  the  surf  zone.  The  surf  zone  is  a  transient  entity  which  varies  in  width  off- 
shore depending  on  the  local  sea  conditions.  These  local  conditions  include  the  beach 
profile,  sea  severity  and  direction  of  wave  field.  In  order  to  define  the  surf  zone  in  this 
study,  wave  field  data  were  analyzed  which  represented  the  range  of  sea  conditions 
experienced  over  the  two  day  period.  The  first  condition,  representing  the  mildest  sea 
condition,  was  taken  as  2155,  23  October  1980.  The  second  condition,  representing  the 
most  severe  sea  condition  was  taken  as  1315,  25  October  1980. 

Linear  wave  theory  will  be  used  to  determine  the  effects  of  shoaling  as  the  wave  field 
enters  shallow  water.  Refraction  will  not  be  considered  since  this  process  tends  to  reduce 

the  shoaling  wave  heights  and  thus  move  the  surf  zone  shoreward.  This  is  shown  by  the 

|cos#0 

definition  of  the  refraction  coefficient,  Kr  =  J where  0O  and  0  are  the  deep  water 

V  cos# 

wave  angle  and  the  finite  water  depth  wave  angle  respectively.  This  ratio  will  always  be 

less  than  one  causing  a  reduction  in  the  resulting  wave  height.  In  order  to  ensure  that 

conservative  results  are  obtained,  refraction  will  be  neglected  in  this  study.  The  beach 

profile  at  the  ARSLOE  project  site  shows  slightly  varying  beach  slopes  so  a  wave 

breaking  criteria  which  considers  the  beach  slope  will  be  used.  Weggel  (1972) 

reinterpreted  laboratory  results  and  showed  a  dependency  of  breaker  height  on  beach 

slope.  The  following  equations  were  used  to  define  the  surf  zone  locations: 


Conservation  of  Energy:  H  2    =    H  x 


C., 
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Breaker  Criterion:     Hb  =  Kdb   ,  where  k  is  defined  by  Weggel  as  follows: 
K  =  b{m)-a(m)^\  with      a(m)  =  43.8(L0-e"19m) 

b(m)  =  156(10 +  e~l9imyl 
As  the  beach  slope,  m  approaches  zero,  this  equation  reduces  to  k  =  0.78. 


The  Group  Velocity,  Cg  =  [  —  I  — 


f  2kd 

1  h 1     will  be  determined  based  on  the 

^      sinh(2Atf), 

dispersion  relation,  I  —     =  gk  tanh(kd) ,  where  L  =  wave  length,  T  =  wave  period 

(assumed  constant),  g  =  gravitational  acceleration,  k  =  the  wave  number  and  d  =  local 
water  depth. 

In  order  to  standardize  the  location  of  the  surf  zone,  the  root-mean- square  wave 
height,  Hrms ,  will  be  used  to  represent  the  irregular  wave  field.    H  i 


=M"' 


where  H;  is  defined  as  the  wave  height  of  the  i   wave  of  the  irregular  wave  field.  The  surf 
zone  which  results  from  the  use  of  H,™  will  represent  the  approximate  "center  of  mass" 
location  where  the  most  significant  set  of  individual  waves  in  the  wave  field  begin  to 
break.  The  surf  zone  calculations  were  based  on  a  root-mean- square  wave  height  that  was 
determined  from  the  deep  water  wave  displacement  data.  This  single  wave  height  was 
used  to  represent  the  irregular  sea  and  shoaled  landward.  Since  the  real  sea  will  have  a 
significant  number  of  waves  with  heights  larger  and  smaller  than  the  root-mean-square,  the 
actual  breaking  zone  is  an  obscure  entity.  This  technique  for  determining  the  surf  zone 
locations  is  adequate  to  make  general  comparisons  when  analyzing  the  spatial  changes  in 
the  nonlineariry  of  the  irregular  sea  in  finite  water  depths.  The  surf  zone  locations  during 
the  mildest  and  most  severe  sea  conditions  over  the  two  day  period  are  shown  in  Table 
3.2. 

In  Table  3.2  it  is  shown  that  during  the  mild  sea  condition,  the  breakpoint,  indicating 
the  beginning  of  the  surf  zone,  is  located  70  meters  from  the  shoreline  in  1.5 1  meters  of 
water.  This  breakpoint  is  the  location  where  a  wave  having  a  height  equivalent  to  the 
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root-mean- square  wave  height  of  the  wave  field  will  break.  During  the  severe  sea 
condition,  this  breakpoint  occurs  1 10  meters  from  the  shoreline  in  3.32  meters  of  water. 
Realistically,  waves  would  be  observed  breaking  considerably  seaward  of  the  breakpoints 
determined  in  the  respective  sea  conditions  due  to  the  spectrum  of  wave  heights  present  in 
an  irregular  sea.  The  root-mean-square  wave  height  breakpoint  is  used  in  order  to 
standardize  the  position  of  the  surf  zone  for  camparison  purposes. 

Table  3.2  Surf  Zone  Locations  for  Mild  and  Severe  Sea  Conditions 


Sea  Condition 

Depth  where  a  wave  with 
height  Hrms,  Breaks 

Start  of  Surf  Zone 
-Dist.  from  Shore- 

Mild  -  10/23  @  2155 
Severe  -  10/25  @  1315 

1.51m 
3.32  m 

70  m 
110m 

3.3    Non-Gaussian  Probability  Density  Function 

The  probabilistic  approach  used  in  this  study  will  make  use  of  the  probability 
density  function  applicable  to  non-gaussian  random  processes  developed  by  Ochi  and  Ahn 
(1994a).  This  probability  density  function  appears  to  represent  non-Gaussian  waves  in 
shallow  and  intermediate  water  depths  as  well  as  Gaussian  waves  in  deep  water. 

3.3.1    Non-Gaussian  Probability  Density  Function  Background 

The  probability  density  function  has  a  background  of  the  Kac-Siegert  solution 
developed  for  the  response  of  non-linear  systems.  The  density  function,  y(t),  is  given  as 
the  sum  of  the  standardized  normal  variates  and  its  squares  as  follows 
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The  parameters  |3j  and  X,j  are  determined  by  finding  the  eigenfunction  and  eigenvalues 

of  the  integral  equation  given  by 

J  K(a>19a>2)y  j(o)2)da)2  =  X  ^  j(eo  J 

where  K(coi  ,02)  =  H(coi  ,  o>2  )  yjSio^JSio^) 

\|/j(co)  =  orthogonal  eigenfunction  =  {    j  y/ j(co)y/ k(k)dco  =  1  for  j  =  k 

otherwise,  vj/j(co)  =  0 
H(coi  ,0)2)  =  second  order  frequency  response  function 
S(co)  =  output  spectral  density  function 
As  can  be  seen  from  the  form  of  this  solution,  the  probability  density  function  of  y(t), 
can  not  be  found  without  knowledge  of  the  second  order  frequency  response  function, 
H(coi  ,  o>2 )  and  there  is  no  way  to  evaluate  this  response  function  for  irregular  seas.  Ochi 
and  Ann  showed  however,  that  Kac-Siegerts'  solution  can  be  obtained  from  the  time 
history  of  the  wave  record  without  any  knowledge  of  the  second  order  frequency  response 
function,  H(coi  ,  oh  ).  The  probability  density  function  thus  obtained  for  describing  non- 


gaussian  random  processes  is 

1 


f(x)  =      ; — exp 

\2nc7 . 


(1  -  ya/jm  -  exp(-yax))2  -  yax 


2(yaamy 

,   L28/cv-:x>0 
where,  y  =  { 

1  3.00/or:x<0 

The  parameters,  a,  fi*  and  o2*  are  determined  from  the  cumulants  evaluated  from  the  wave 
data.  These  cumulants  are  related  to  the  Pj  and  X}  values  of  the  Kac-Siegert  solution  by 
the  following  relationships 

2  N 


2  N  2  N 

k,  =  6^  J32JAJ  +sfdXij  ,etc. 
7=1  j=i 
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In  the  present  study,  ki  represents  the  mean  water  level  and  is  therefore  defined  by, 
ki  =  0.  k2  and  k3  are  equal  to  the  data  sample  moments,  E[x2]  and  E[x3]  respectively. 
Based  on  this  relationship,  the  parameters,  a,  u*  and  a2*  can  be  determined  through  the 
iterative  solution  of  three  equations  along  with  the  analysis  of  a  data  record  for  the  non- 
Gaussian  random  process  being  analyzed.  The  three  equations  are 

a  a2*  +  a  |i*2  +  u*  =  E[x]  =  0 

a2.-2aV*  =  E[x2] 

2a  a4*(3  -  8aV.)  =  E[x3] 

3.3.2    Probability  Density  Function  Application  to  Wave  Data 

In  order  to  determine  values  for  the  three  parameters,  a,  p.*  and  a"* ,  which  are 
used  to  define  the  degree  of  non-linearity  a  shoaling  wave  field  has  at  a  given  location  and 
sea  state,  data  records  must  be  analyzed.  With  the  applicability  of  the  probability  density 
function  verified  based  on  the  wave  data  collected  during  ARSLOE,  the  approximately 
1000  data  records  from  1200  on  23  October  to  13 15  on  25  October  will  be  analyzed  to 
determine  values  for  the  three  parameters  defining  the  probability  density  function.  The 
procedure  used  during  this  study  to  perform  this  task  is  as  follows: 

-  Compute  the  standard  deviation  of  each  twenty  minute  wave  record  consisting  of 
4800  discrete  wave  displacement  data  points  E(x2 )  =  —  V  x2,-  . 

(E[x]  does  not  need  to  be  computed  since  it  is  defined  as  zero  and  the  data  are 
already  standardized). 


-  Compute  the  third  moment  of  the  same  twenty  minute  wave  record 
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-  Define  the  standard  deviation  of  the  data  record  as  o  =  y]E(x2 )  . 

-  Evaluate  the  three  parameters,  a,  u»  and  a2* ,  through  iterative  solution  of  the 
following  three  equations: 

a  ct2.  +  a  u»2  +  |x*  =  E[x]  =  0 

a2.-2aV.  =  E[x2] 

2a  ct4*(3  -  8aV.)  =  E[x3]. 
It  is  important  to  note  that  the  standard  deviation  of  the  data  record  is  used  in  the 
analysis  as  an  indicator  of  the  sea  severity  and  not  the  significant  wave  height.  This  is 
because  the  definition  of  significant  wave  height  is  not  applicable  to  non-Gaussian  waves. 

All  available  data  sets  covering  the  seven  finite  water  depth  gage  locations 
between  1200,  23  October  1980  and  1400,  25  October  1980  were  analyzed  and  a,  u*  and 
g2»  were  determined. 

Based  on  a  plot  of  the  individual  parameters  in  the  space  domain,  it  was 
determined  that  a  smoothing  technique  should  be  used  to  refine  the  "a"  parameter  values. 
This  smoothing  was  done  in  order  to  remove  the  sporadic  tertiary  effects  while  ensuring 
the  smoothed  "a"  value  does  not  lose  any  authenticity.  An  exponential  smoothing 
technique  was  used  which  resulted  in  the  lowest  mean  absolute  deviation  (MAD)  between 
the  raw  "a"  parameter  value  and  the  smoothed  "a"  parameter.  The  remaining  two 
parameters,  u*  and  o  *  did  not  show  the  same  susceptibility  to  tertiary  effects  so  no 
smoothing  was  done.  Comparisons  between  the  actual  and  smoothed  "a"  parameters  are 
shown  in  Appendix  A  for  all  Gage  positions. 

3.3.3    Dimensional  Analysis  of  PDF  Parameters 

For  each  of  the  approximately  one-hundred-fifty  data  sets  per  gage  position,  the 
probability  density  function  parameters,  a,  \i* ,  a*  were  determined  iteratively.       Each  of 
the  parameters  contain  discrete  dimensions.  Based  on  dimensional  analysis,  these 
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dimensions  are  as  follows:  It  is  known  that  k2,  representing  the  data  sets'  variance  has 
dimensions  of  [L2].  It  is  also  known  that  k3,  representing  the  data  sets'  third  moment  has 
dimensions  of  [L3].  It  follows  from  the  knowledge  of  these  dimensions  that  the 
parameters,  a,  u* ,  a*  must  have  dimensions  [1/L],  [L]  and  [L]  respectively. 

In  order  to  minimize  any  dependence  on  external  conditions  and  to  ensure  the 
solutions  are  universally  applicable,  these  parameters,  a,  \i* ,  a*  must  be  made  non- 
dimensional.  The  non-dimensional  quantities  that  will  be  used  in  this  study  are,  aa,  u*/a, 
a*/a  and  a/d.  a  is  defined  as  the  standard  deviation  of  that  wave  displacement  data  and  d 
is  defined  as  the  water  depth. 


CHAPTER  4 
PROBABILITY  DENSITY  FUNCTION  ANALYSIS 


4.1    Probability  Density  Function  Applicability  Verification 

As  stated  in  the  introduction,  one  of  the  objectives  of  this  study  is  to  determine 
whether  or  not  the  probability  density  function  developed  by  Ochi  and  Ahn  is  valid  for 
representing  waves  in  finite  water  depths.  To  verify  the  applicability  of  this  probability 
density  function,  wave  displacement  data  sets  obtained  from  the  ARSLOE  experiment 
were  used  in  comparing  the  histogram  of  the  data  record  to  the  probability  density 
function.  Each  data  set  contained  4800  wave  displacement  readings  taken  over  a  twenty 
minute  period.  To  ensure  a  thorough  verification  of  the  probability  density  function,  the 
data  records  used  for  the  analysis  were  from  all  seven  "finite"  water  depth  gage  locations 
and  analysis  was  done  over  a  wide  range  of  sea  state  conditions.  One  analysis  was  done 
on  data  from  23  October  1980  at  2155.  This  represented  a  mild  sea  state.  The  other 
analysis  was  done  on  data  from  25  October  1980  at  1315.  This  data  set  was  at  the  peak 
of  the  most  severe  storm  experienced  during  the  entire  ARSLOE  study. 

Figures  4. 1-4. 14  show  comparisons  between  the  histogram  data  and  the  non- 
Gaussian  probability  density  function.  Included  also  in  the  figures  is  the  Gaussian 
(normal)  probability  density  function  for  comparison.  Figures  4. 1-4.7  pertain  to  the  mild 
sea  condition  while  Figures  4.8-4. 14  represent  the  severe  sea  condition.  As  can  be  seen 
from  these  figures,  the  probability  density  function  agrees  very  well  with  the  histogram 
generated  from  the  data  record  in  all  cases.  This  implies  that  the  probability  density 
function  derived  by  Ochi  and  Ahn  can  be  used  with  confidence  when  analyzing  non- 
Gaussian  random  processes  such  as  wave  displacements  in  shallow  and  intermediate 
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waters.  It  can  also  be  seen  that  the  probability  density  function  approximates  the  normal 
distribution  when  the  wave  field  is  in  the  deeper  intermediate  water  condition  under  mild 
sea  conditions. 
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Table  4.1  Gage  615  -  10/23  @  2155 


WAVE  PROPERTIES 


615 

10/23  (a),  2155 

a  (1/m) 

.572645 

Mm) 

- .04948 

a.  (m) 

.275683 

<T(m) 

.267219 

d(m) 

1.74 

CT/d 

.153574 

Figure  4.1  Histogram  Data 
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Table  4.2  Gage  635  -  10/23  @  2155 


WAVE  PROPERTIES 


635 

10/23  (a).  2155 

a  (1/m) 

.481441 

Mm) 

-.0547 

a.  (m) 

.318871 

CT(m) 

.30989 

d(m) 

2.83 

a/d 

.109502 

Figure  4.2  Histogram  Data 
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Table  4.3  Gage  645  -  10/23  @  2155 


WAVE  PROPERTIES 


645 

10/23  dab.  2155 

a  (1/m) 

.200988 

M.  (m) 

-  .03084 

a.  (ra) 

.343446 

a(m) 

.340708 

d(m) 

4.31 

a/d 

.079051 

Figure  4.3  Histogram  Data 
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Table  4.4  Gage  655  -  10/23  @  2155 


WAVE  PROPERTIES 


655 

10/23  (a)  2155 

a  (1/m) 

.104437 

Mm) 

- .01771 

a.  (m) 

.340056 

a(m) 

.339137 

d(m) 

6.0 

a/d 

.056523 

Table  4.5  Gage  665  -  10/23  @  2155 


665 
a  (1/m) 
H.  (m) 
ct.  (m) 
a(m) 
d(m) 
a/d 


WAVE  PROPERTIES 
10/23  deb,  2155 

.126266 

-  .01835 
.365155 
.364236 

6.08 
.059907 


Figure  4.5  Histogram  Data 
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Table  4.6  Gage  675  -  10/23  @  2155 


WAVE  PROPERTIES 


675 

10/23  (2)  2155 

a  (1/m) 

.094349 

Mm) 

-  .01014 

ct.  (m) 

.375221 

fj(m) 

.374948 

d(m) 

6.82 

cr/d 
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Figure  4.6  Histogram  Data 
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Table  4.7  Gage  625  -  10/23  @  2155 


625 
a  (1/m) 
M.  (m) 
a.  (m) 
o-(m) 
d(m) 

a/A 


WAVE  PROPERTIES 

10/23  (a)  2155 
.069914 

-  .00819 
.394571 
.394401 

9.66 
.040828 


Figure  4.7  Histogram  Data 
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Table  4.8  Gage  615  -  10/25  @  1315 


WAVE  PROPERTIES 


615 

10/25  (a),  1315 

a  (1/m) 

.372027 

Mm) 

-  .09863 

a.  (m) 

.498230 

cr(ni) 

.47981 

d(m) 

2.23 

a/A 
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Figure  4.8  Histogram  Data 
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Table  4.9  Gage  635  -  10/25  @  1315 


WAVE  PROPERTIES 


635 

10/25  (2)  1315 

a  (1/m) 

.318011 

H.  (m) 

- .10497 

ct,  (m) 

.591069 

a(m) 

.573397 

d(m) 

3.32 

a/d 

.17271 

Figure  4.9  Histogram  Data 
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Figure  4.10  Histogram  Data 


Table  4.10  Gage  645-10/25  @  1315 


WAVE  PROPERTIES 


645 

10/25  (a).  1315 

a  (1/m) 

.256391 

H.  (m) 

-  .10245 

a.  (m) 

.628243 

CT(m) 

.612186 

d(m) 

4.66 

a/d 
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Table  4.11  Gage  655  -10/25  @  1315 


WAVE  PROPERTIES 


655 

10/25  (a)  1315 

a  (1/m) 

.269927 

M.  (m) 

-  .12356 

a.  (m) 

.658813 

a(m) 

.636821 

d(m) 

6  .44 

a/d 

.098885 

Figure  4.11  Histogram  Data 
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Table  4.12  Gage  665  -10/25  @  1315 


WAVE  PROPERTIES 


665 

10/25  (2)  1315 

a  (1/m) 

.269743 

Mm) 

-  .12978 

a.  (m) 

.67434 

a(m) 

.650698 

d(m) 

6.52 

a/d 

.09980 

Figure  4.12  Histogram  Data 
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Table  4.13  Gage  675  -10/25  @  1315 


WAVE  PROPERTIES 


675 

10/25  (2)  1315 

a  (1/m) 

.220143 

H.  (m) 

-  .12062 

<r.  (m) 

.727669 

<r(m) 

.708477 

d(m) 

7.22 

a/d 

.098127 

Figure  4.13  Histogram  Data 
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Table  4.14  Gage  625  -10/25  @  1315 


WAVE  PROPERTIES 


625 

10/25  (a).  1315 

a  (1/m) 

.140009 

H.  (m) 

-.10961 

tr.  (m) 

.872244 

<j(m) 

.858792 

d(m) 

10.04 

<r/d 

.085537 

Figure  4.14  Histogram  Data 


Based  on  a  comparison  with  the  water  depth  regions  defined  in  Section  3.2. 1,  it 
can  be  seen  that  the  histogram  data  and  the  probability  density  function  closely  resemble 
the  normal  distribution  even  after  the  wave  field  enters  intermediate  water  for  milder  sea 
conditions.  This  is  shown  in  Figures  4.6  and  4.7.  This  implies  that  the  characteristics  of 
the  wave  profile  are  dependent  on  the  sea  condition  as  well  as  the  water  depth. 

The  histogram  data  for  both  the  mild  and  severe  sea  conditions  in  Figures  4. 1-4. 14 
show  the  distribution  of  the  wave  profiles  progressively  becomes  more  non-Gaussian  as 
the  waves  propagate  toward  shallower  depths.  The  non-Gaussian  trend  becomes  more 
extreme  as  the  sea  condition  becomes  more  severe.  It  is  significant  to  note  also  that  the 
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non-Gaussian  characteristics  are  most  pronounced  at  the  shallower  gage  locations,  Gages 
615  and  635. 

4.2     Gaussian/non-Gaussian  Boundary  Definition 

The  histogram  data  shown  in  Figures  4. 1  -  4. 14  represent  the  same  wave  profiles 
depicted  in  Figures  3.2  -  3. 15.  The  transformation  from  Gaussian  to  non-Gaussian 
profiles  begin  with  subtle  changes  in  the  amplitude  of  the  negative  displacements  and  the 
width  of  the  troughs.  This  initial  transformation  can  be  seen  be  comparing  Figures  4.4  and 
4.6.  Figure  4.6  shows  a  nearly  normal  distribution.  Figure  4.4  shows  a  slight  elevation  of 
the  apex  and  a  narrowing  of  the  negative  tail.  It  is  also  significant  to  note  that  the  positive 
part  of  the  distribution  remains  relatively  unchanged.  This  indicates  that  the  initial  changes 
affect  the  trough  region  of  the  profile  more  drastically  than  the  positive  region.  Figure  4.4 
demonstrates  that  the  trough  region  becomes  shallower  and  wider.  As  the  profile 
continues  to  shoal,  changes  begin  to  occur  in  the  positive  displacement  region  of  the 
distribution.  The  peaks  become  narrower  and  higher.  This  can  be  seen  graphically  by 
comparing  Figures  4.2  and  4.4.  The  positive  region  of  the  wave  profile  is  not  as  full  when 
compared  with  the  normal  distribution  and  the  positive  tail  extends  farther  than  the  normal 
curve. 

As  described  above,  significant  transformations  occur  in  the  wave  profile  as  the 
wave  train  propagates  into  shallower  water.  There  is  a  point  in  this  transformation  when 
the  wave  profile  deviates  significantly  enough  from  the  normal  distribution  that  the  waves 
can  no  longer  be  considered  a  Gaussian  random  process.  Identification  of  this  boundary  is 
important  to  ensure  assumptions  made  regarding  the  waves  remain  valid.  Analysis  of  the 
histogram  data  shows  that  the  non-Gaussian  characteristics  of  coastal  waves  depends  on 
water  depth  and  sea  severity.  Based  initially  on  analysis  of  the  histogram  data  shown  in 
Figures  4.5  and  4.6,  it  appears  this  transition  boundary  occurs  when  the  non-dimensional, 
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wave  standard  deviation,  a/d,  is  approximately  0.06.  The  non-dimensional  parameter, 
a/d,  is  used  because  it  considers  both  sea  severity  and  water  depth.  Eight  random  sea 
conditions  were  selected  having  a/d  values  between  0.04  and  0. 10.  These  eight  data  sets 
are  depicted  in  Figure  4. 15  in  terms  of  aa  as  a  function  of  a/d  .  The  data  points  labeled  A- 
H  in  Figure  4.23,  each  correspond  to  a  distribution  shown  in  Figures  4. 15-4.23.  Data 
points  A-G  were  intentionally  grouped  around  the  10  meter  water  depth.  This  minimized 
any  impact  water  depth  and/or  beach  slope  may  have  on  the  analysis  results.  The 
statistical  data  for  these  points  are  provided  in  Table  4. 15. 


Figure  4.15  Sample  Data  Points  to  Analyze  Gaussian/Non-Gaussian  Boundary 
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Table  4.15  Statistical  Properties  of  Data  Points  A-H  of  Figure  4.15 


POINT 

a/d 

H, 

Water  Depth,  d 

A 

0.1000 

3.08  m 

9.51m 

B 
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2.74  m 

9.67  m 

C 
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D 
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Figure  4.16  Point  A  of  Figure  4.15 
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Figure  4.17  Point  B  of  Figure  4.15 
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Figure  4.18  Point  C  of  Figure  4.15 
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Figure  4.19  Point  D  of  Figure  4.15 
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Figure  4.21  Point  F  of  Figure  4.15 
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Figure  4.22  Point  G  of  Figure  4.15 
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Figure  4.23  Point  H  of  Figure  4.15 


Figures  4. 16  -  4.23  confirm  that  the  region  between  o7d  =  0.06  and  o7d  =  0.08, 
shown  in  Figures  4.20  and  4. 17  respectively,  appears  to  be  when  the  wave  profile 
transforms  from  Gaussian  to  non-Gaussian.  A  closer  inspection  of  this  region  is  shown  in 
Figures  4. 18  -  4.20.  The  distribution  shows  a  considerable  increase  in  deviation  from  the 
normal  between  o7d  =  0.06  and  o7d  =  0.0616  as  shown  in  Figures  4.20  and  4. 19 
respectively. 

Due  to  the  subjective  nature  of  this  boundary,  a  general  boundary  can  be  defined 
based  on  probabilistic  analysis.  This  boundary,  defining  the  conditions  under  which  the 
wave  profile  can  no  longer  be  considered  Gaussian,  occurs  when  o7d  =0.06. 

A  similar  study  was  conducted  by  Ochi  and  Wang  (1984);  however,  they  defined 
the  boundary  condition  between  Gaussian  and  non-Gaussian  wave  profiles  in  terms  of  the 
wave  record  parameter,  ^3.   X3  represents  the  skewness  of  the  wave  profile  which  they 
determined  to  be  the  dominant  parameter  affecting  the  non-Gaussian  characteristics  of 
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coastal  waves.  Ochi  and  Wang  concluded  that  coastal  waves  in  seas  for  which  X3  is  less 
than  0.2  can  be  considered  as  a  Gaussian  random  process.  For  comparison  purposes,  this 
boundary  condition  is  plotted  along  with  the  findings  of  this  present  study  in  Figure  4.24. 
In  terms  compatible  with  the  approach  of  this  present  study,  their  boundary  condition  can 
be  defined  as  a  /d  =  0.0543.  It  is  encouraging  that  there  shows  very  close  agreement 
between  both  findings.  Additional  analysis  was  contributed  by  Ochi  and  Ann  ( 1 994b) 
whose  data  points  are  shown  in  Figure  4.24  as  white  circles.  Ochi  and  Ahn  (1994b) 
defined  the  transition  criterion  between  Gaussian  and  non-Gaussian  as  occurring  when  a  < 
0.01  where  "a"  is  one  of  the  parameters  in  the  non-Gaussian  probability  density  function. 
The  white  circle  data  points  define  the  limiting  cases  based  on  his  criterion.  This  third 
method  for  determining  the  boundary  between  Gaussian  and  non-Gaussian  wave  profiles 
agrees  very  well  with  the  other  two  further  verifying  the  validity  of  this  finding. 
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CHAPTER  5 
TREND  ANALYSIS  OF  PARAMETERS 


5.1    Relationship  between  a  and  Water  Depth 

As  was  discussed  earlier,  the  distribution  of  the  wave  profile  is  highly  dependent 
on  sea  severity  and  water  depth  with  sea  severity  being  defined  in  terms  of  the  wave 
profiles'  standard  deviation,  a.  Since  these  two  terms  will  be  used  considerably  in  the 
following  trend  analyses,  their  relationship  with  each  other  will  be  examined. 

Figure  5.1  shows  examples  of  simultaneous  measurements  of  a  plotted  as  a 
function  of  distance  offshore  for  two  cases.  It  is  evident  from  Figure  5.1,  that  the 
standard  deviation  of  the  wave  profile  is  influenced  by  water  depth.  As  water  depth 
decreases,  the  relative  magnitude  of  the  standard  deviation  decreases.  It  is  also  seen  that 
for  more  severe  sea  conditions,  represented  by  10/25  @  1315  data  in  Figure  5. 1,  the 
overall  change  in  ct  as  the  wave  field  propagates  into  shallower  waters  is  greater  as 
compared  with  milder  sea  conditions  (10/23  @  2155).  Based  on  theoretical 
considerations  of  o,  it  can  be  assumed  that  a  approaches  zero  as  distance  offshore 
approaches  zero.  This  trend  is  observed  in  Figure  5. 1  since  both  cases  shown  appear  to 
converge  to  the  origin.  From  the  information  on  the  surf  zone  locations  for  these  two 
cases  computed  in  Chapter  3,  it  appears  that  o  begins  to  decrease  rapidly  towards  zero 
once  the  break  point  is  reached. 

The  information  shown  in  Figure  5. 1  based  on  two  cases  is  shown  in  Figure  5.2  for 
all  data  sets  analyzed  in  this  study.  Instead  of  a  being  plotted  as  a  function  of  distance 
offshore,  Figure  5.2  plots  a  as  a  function  of  water  depth.  This  will  take  into  consideration 
any  local  changes  in  mean  sea  level. 
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Plot  of  o  vs  Distance  Offshore 
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Figure  5.1  Spatial  Comparison  of  a  at  Mild  and  Severe  Sea  Conditions 

The  trends  evident  in  Figure  5.2  reinforce  the  analysis  done  for  the  two  cases  in 
Figure  5.1.  The  upper  bounds  of  the  discrete  data  represent  the  severe  sea  conditions 
while  the  lower  bounds  represents  the  mild  sea  conditions. 
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Figure  5.2  All  a  Data  Plotted  as  a  Function  of  Water  Depth 
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It  can  also  be  inferred  from  Figure  5. 1  and  Figure  5.2  that  as  the  sea  state  becomes 
milder,  the  waves'  standard  deviation  approaches  a  constant  value  and  hence  the  sea 
severity  approaches  a  constant  before  decreasing  towards  the  origin. 

In  the  absence  of  elaborate  wave  data,  the  value  for  a  can  be  readily  determined 


with  a  knowledge  of  the  local  and  breaking  wave  heights  and  the  relative  water  depth, 
d_ 

gT 


defined  by  — —  ,  where  g  =  gravitational  constant  and  T  =  wave  period,  using  a  method 


such  as  stream  function  theory  (Dean  1974).  To  ensure  the  use  of  accurate  wave  height 
values  in  applying  the  stream  function  theory,  nonlinear  shoaling  methods  should  be  used 
in  finite  water  depths  to  determine  the  breaking  and  local  wave  heights. 

5.2    Non- Gaussian  PDF  Parameter  Analysis 

The  non- Gaussian  probability  density  function  parameters,  a,  \i*  and  a,,  will  be 
analyzed  in  non-dimensional  form  spatially  as  a  function  of  distance  offshore,  against  each 
other  in  non-dimensional  form,  and  with  aa  as  a  function  of  o/d.  The  first  analysis  of  the 
non-Gaussian  probability  density  function  parameters  will  involve  comparing  the  non- 
dimensional  parameters  against  the  distance  the  gage  was  located  offshore.  This  will 
provide  a  spatial  orientation  of  the  parameters  which  will  disregard  any  local  fluctuations 
in  water  depth  due  to  tidal  cycle  of  storm  surge.  The  results  are  shown  in  Figures  5.3-5.5. 

5.2.1    Spatial  Analysis  of  Parameters 

All  three  parameters  shown  in  Figures  5.3-5.5  appear  to  gradually  increase  in 
absolute  magnitude  until  approximately  150  meters  to  175  meters  offshore.  The  mild 
condition  represented  by  10/23  @  2155  consistently  yields  lower  values  when  compared 
with  the  more  severe  sea  condition  on  10/25  @  1315.  As  the  sea  condition  becomes  more 
mild,  the  wave  profile  resembles  the  normal  distribution.  When  the  wave  profile 
represents  exactly  the  normal  distribution,  "a"  equals  zero  and  a*  equals  a.  It  can  be  seen 
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from  Figures  5.3  and  5.5  that  as  approaches  zero  under  mild  conditions  and  s*/s 
approaches  one.  Figure  5.4  shows  that  m»/s  approaches  zero  under  these  same  conditions 
which  is  also  consistent  with  the  theoretical  limit  of  the  non-Gaussian  probability  density 
function. 

There  is  also  a  significant  perturbation  evident  in  Figures  5.3-5.5  around  150 
meters  offshore.  The  mild  sea  condition  shows  this  abrupt  change  occurring  seaward 
when  compared  with  its  location  under  severe  sea  conditions  and  the  perturbation  is  more 
pronounced  in  the  severe  sea  condition  case.  In  general,  all  trends  analyzed  for  the  non- 
dimensional  parameter,  m*/s,  are  almost  identical  results  to  those  found  for  as  and  s»/s  but 
in  the  mirror  image.  This  is  due  to  the  fact  that  the  value  of  m*  is  negative.  The  negative 
sign  is  supported  by  the  theory  used  to  develop  the  probability  density  functions'  three 
governing  equations. 

The  next  spatial  analysis  will  consist  of  comparing  the  non-dimensional  parameters 
with  water  depth.  This  comparison  will  be  sensitive  to  any  influence  of  local  depth 
changes  occurring  from  tidal  cycles  and/or  storm  surges.  Figures  5.6  -  5.8  show  the 
results  of  this  analysis. 

The  trends  from  this  analysis  are  the  same  as  those  determined  by  Figures  5.3-5.5. 
Due  to  the  difference  in  scale  of  the  abscissa  the  curves  based  on  water  depth  appear  to 
show  much  smoother  transitions.  Figures  5.6-5.8  include  two  additional  cases 
representing  intermediate  sea  conditions  between  those  represented  by  10/23  @  2155  and 
10/25  @  1315.  These  two  additional  cases  approximate  the  trends  displayed  by  the  mild 
sea  condition  of  10/23  @  2155.  This  is  consistent  with  the  data  and  time  of  these  cases 
since  they  occurred  prior  to  the  severe  storm  so  the  sea  condition  was  not  yet  impacted 
significantly. 

As  was  discussed  earlier,  in  the  deeper  water  conditions  all  three  parameters  tend 
to  show  convergence  towards  their  respective  values  for  a  Gaussian  wave  profile.  In 
extremely  shallow  water  conditions,  the  parameters  also  tend  to  converge  towards  a 
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constant  value  which  is  different  depending  on  the  sea  condition.  For  the  milder  sea 
conditions,  the  parameters  tend  to  converge  to  the  values  shown  in  Table  5.1. 


Table  5.1  General  Shallow  Water  Limits  of  Parameters 


Non-Dimensional  Parameter 

Mild  -  10/23  @  2155 

Severe  -  10/25  @  1315 

aa 

0.15 

0.18 

u./a 

-0.18 

-0.22  to  -0.25 

oJg 

1.033 

1.04  to  1.045 

When  the  separate  plots  for  each  non-dimensional  parameter  are  compared  in  both 
spatial  cases,  it  is  apparent  that  similar  trend  descriptions  fit  all  three  parameters.  For  this 
reason,  a  closer  analysis  will  be  made  among  the  three  variables  to  determine  if  this 
similarity  can  be  quantified.  The  next  trend  analysis  consists  of  comparing  the  three 
parameters  with  each  other. 


5.2.2    Relationship  Between  the  three  non-dimensional  parameters,  aa,  u»/q  and  aja 

As  discussed  earlier,  the  three  non-dimensional  parameters  appear  to  be  highly 
correlated  with  each  other.  A  closer  look  will  be  taken  to  see  if  specific  relations  can  be 
developed  connecting  the  three  non-dimensional  parameters  together.  Figure  5.9  and 
Figure  5. 10  show  plots  of  the  raw  data  representing  u»/o  and  cr./a  compared  with  acr 
respectively  at  the  extreme  water  depth  locations,  gage  615  and  gage  625. 

All  data  sets  from  Gages  615  and  625  were  plotted  since  they  represent  the 
shallowest  and  deepest  gage  locations.  There  are  very  definite  relationships  which  hold 
constant  between  the  extreme  range  of  water  depths  and  sea  conditions,  that  link  ao  to 
\x.J<3  and  cJg.  The  functional  relationship  can  be  determined  only  after  considering  all 
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the  data  points  from  all  water  depth  and  sea  severity  conditions  to  ensure  it  remains 
consistently  valid.  A  deterministic  relationship  between  these  non-dimensional  parameters 
is  a  topic  of  active  research. 
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Figure  5.9  Plot  of  (li*/ct  as  a  function  of  aa  for  Gages  615  and  625 
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Figure  5.10  Plot  of  a»/a  as  a  function  of  aa  for  Gages  615  and  625 
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5.2.3    Relationship  Between  aa  and  a/d 

Based  on  the  apparent  correlation  between  the  parameters  of  the  non-Gaussian 
probability  density  function,  aa  will  be  used  to  represent  the  general  trends  of  the  other 
two  parameters  when  plotted  as  a  function  of  a/d.  This  relationship  is  provided  in  Figure 
5. 1 1  to  demonstrate  general  trends  between  the  wave  aa,  and  the  local  environment 
represented  in  terms  of  sea  severity/water  depth  in  non-dimensional  form  The  lower 
values  of  a/d  correspond  to  the  deeper  gage  locations  and  the  trend  is  for  these  values  to 
approach  the  values  defined  by  wave  profiles  having  a  normal  distribution.  This  plot  also 

t 

reinforces  the  analysis  done  in  Chapter  4  which  determined  the  boundary  values  of  the 
Gaussian  distribution.  At  a/d  =  0.06,  the  aa  values  appear  to  begin  converging  towards 
the  Gaussian  limit  located  at  the  origin  of  Figure  5.11. 

The  higher  values  of  a/d  represent  the  shallower  water  depths.  It  appears  that  at 
these  values  aa  tends  to  approach  a  finite  asymptotic  limit.  This  is  consistent  with  the 
findings  from  the  spatial  analysis  figures.  In  between  these  two  limits  of  Figure  5.11  there 
is  a  depression  in  the  severe  sea  condition  case.  This  phenomenon  must  be  analyzed  more 
thoroughly  to  ascertain  a  reasonable  explanation  for  its  existence  and  will  be  the  subject  of 
a  later  study. 
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CHAPTER  6 
SUMMARY  AND  CONCLUSIONS 


It  is  shown,  based  on  the  results  of  this  study,  that  wind  generated  waves  exhibit 
significant  changes  in  their  profile  forms  as  they  propagate  into  finite  water  depths.  The 
actual  nature  of  this  transformation  has  been  a  subject  that  has  not  been  covered  in  any 
significant  detail.  The  unique  wave  data  records  collected  during  the  ARSLOE  Project  of 
1980  enabled  a  detailed  study  on  this  topic  to  be  produced. 

One  objective  of  this  study  was  to  verify  the  applicability  of  the  non-Gaussian 
probability  density  function  developed  by  Ochi  and  Ahn  (1994b)  defined  as 


f(x)  =      i — — exp 


1  - 

— T(l  -  yafi,  -  exp(-yax))    -  yax 


2(yaa.) 

,   USfor.x>0 
where,  y  =  { 

1  3.00/or:x<0 

Through  statistical  analysis  of  the  ARSLOE  data,  it  has  been  confirmed  that  this 
probability  density  function  can  be  used  with  confidence  to  represent  wave  profiles  in 
finite  water  depths.  It  correlates  very  well  with  histogram  data  over  an  extensive  range  of 
sea  conditions.  The  distribution  of  the  wave  profile  in  finite  water  depths  shows  a 
dependence  on  water  depth  and  sea  severity.  In  particular,  as  water  depth  decreases,  the 
distribution  becomes  increasingly  non-Gaussian.  Also,  as  sea  severity  increases  the  same 
trend  occurs. 

With  the  applicability  of  the  probability  density  function  verified,  an  analysis  of  the 
boundary  condition  where  the  wave  profile  can  no  longer  be  considered  Gaussian  was 
conducted.  A  general  relationship  has  been  proven  to  exist  between  the  largest  significant 
wave  height  and  the  local  water  depth  where  the  wave  profile  can  still  be  considered 
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H  IT 

Gaussian.  This  relationship  is  defined  as  — -  =  0.24  or  —  =  0.06 .  Both  forms  of  this 

relationship  are  identical  since  the  significant  wave  height,  Hs  is  defined  as  four  times  the 
standard  deviation,  a,  in  deep  water  conditions.  This  relationship  has  been  verified  using 
three  different  approaches  with  each  yielding  similar  results.  In  shallower  water 
conditions,  this  boundary,  defining  the  transition  from  Gaussian  to  non-Gaussian  wave 
profiles,  appears  to  shows  a  dependence  on  beach  slope.  This  is  a  qualitative  observation 
from  the  present  study  and  additional  research  must  be  done  to  quantify  this  finding.  This 
dependence  on  beach  slope  may  be  similar  to  the  influence  of  beach  slope  on  breaking 
waves.  The  linear  result  defined  by  this  study  is  akin  to  the  earliest  breaking  criterion 
based  on  Solitary  Wave  Theory. 

At  present,  the  only  way  this  probability  density  function  can  be  defined  is  through 
knowledge  of  statistical  properties  of  the  wave  field  of  interest.  In  particular,  the  mean, 
variance  and  third  moment  of  the  wave  profile  must  be  known.  This  requires  that  wave 
profile  data  be  obtained  in  order  to  evaluate  the  parameters  of  the  probability  density 
function. 

Another  objective  of  this  study  was  to  determine  relationships  linking  the 
parameters  of  the  probability  density  function  to  local  field  conditions.  If  the  wave  profile 
statistics  can  be  defined  in  terms  of  local  environmental  conditions  such  as  sea  state  and 
water  depth,  the  probability  density  function  can  be  defined  without  wave  data.  This  study 
has  revealed  qualitative  trends  between  the  parameters  of  the  probability  density  function 
and  local  field  conditions  based  on  the  extensive  data  complied  from  ARSLOE.  From 
these  qualitative  trends,  it  is  evident  that  the  specific  relationships  between  the  parameters 
and  environmental  conditions  are  complex  and  may  require  consideration  of  additional 
parameters  besides  sea  severity,  c,  and  water  depth.  The  pursuit  of  a  functional 
relationship  between  the  parameters  and  the  local  environment  continues  to  be  a  topic  of 
considerable  interest  and  active  research. 


APPENDIX 
COMPARISON  OF  ORIGINAL  AND  SMOOTHED  "a"  VALUES 


It  can  be  seen  that  the  nonlinearity  parameter,  "a"  is  characterized  by  a  significant 
amount  of  local  flucuations  superimposed  on  the  overall  trend  of  the  parmeter.  Based  on 
the  definition  of  "a"  as  a  nonlinearity  parameter  used  to  describe  the  probability  density 
function  of  a  wave  profile,  these  small  flucuations  can  be  smoothed  without  any  noticable 
loss  of  accuracy.  It  is  highly  unlikely  that  the  probability  density  function  of  a  wave  profile 
changes  so  sporadically  over  such  short  periods  of  time  as  indicated  by  the  minor 
variations  of  "a".  There  was  not  any  one  smoothing  technique  that  was  used  in  the 
smoothing  of  "a"  at  all  gage  positions.  Instead,  several  smoothing  techniques  were 
evaluated  and  the  technique  which  produced  the  lowest  mean  absolute  deviation  (MAD) 
was  the  technique  used  for  that  particular  gage  position. 
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